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Abstract 

We consider 2+1-dimensional classical noncommutative scalar field theory. The general ansatz 
for a radially symmetric solution is obtained. Some exact solutions are presented. Their possible 
physical meaning is discussed. The case of the finite 9 is discussed qualitatively and illustrated by 
some numerical results. 



1 Introduction 



Noncommutative field theories became an object of intense study during the last year due to 
their close connection to string theory. Namely, in the large B-field limit the associative string 
field algebra A factors into an algebra that acts on the string center of mass and an algebra 
that describes all other degrees of freedom[l,2,3]. The second algebra seems to be terribly 
complicated while the first one leads to the effective noncommutative theory. Specifically, 
qq , description of the scalar field theory of tachyons on the world- volume of unstable Dp-branes 

in superstring theory is a field of employment for noncommutative theory [4]. 

Of especial interest are solutions of the equation = <p, i.e. projection operators. 

Here * denotes the nonlocal (Moyal) star product that becomes 

f*g{x) = e^" a >" f{x')g{x")\ x , =x „ =x (1) 
when written in rescaled coordinates. If one goes to the momentum space, 



/(*) = J ^f{x)e- ik *, (2) 

then the ^-product takes the form 

f^gip) = J ^f(k)g(p~k)e-i^ k ^. (3) 

After rescaling the coordinates the (static) energy functional is given by 

E[</>] = 1 J cPx^d^t + 6V{<f)). (4) 

Here -^V((j))\^—o — is implied, 6 is the noncommutativity parameter[5]. All the fields are 
multiplied using the ^-product. 

The trick used to solve the equations of motion is to formally identify the noncommu- 
tative quantum space on which the 0-field lives with the classical phase space of a single 
particle and to recast the problem into the operator formalism using the Weyl-Moyal (WM) 
correspondence. The latter is the correspondence between the multiplication (composition) 
of operators acting on a single-particle Hilbert space and the *-product of functions on the 
corresponding phase space. It is defined by the commutativity of the following diagram: 
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f,g — * fg 
f,g — ► f*g 

n denotes a map (more precisely, one-one correspondence) from the algebra of the functions 
to the algebra of corresponding operators using symmetric (Weyl) ordering prescription: 

: f(p,Q) - f(P,Q) = J ^-f(k P >k q y (kpP+kq4) > ( 5 ) 

f(k) being defined by (2). Thus the ^-product [] provides such a deformation of the algebra of 
functions on the phase space (which is restricted to be a symplectic manifold with a constant 
Poisson bi- vector 6^ v ) that Q, becomes an isomorphism [6]. The following relation also holds: 

^J>(x)=Tr n l (6) 

If the Hilbert space 7i is the space of wave functions in coordinate representation, there 
exists the explicit formula to map a kernel of a Hilbert-Schmidt operator acting on TL into 
the corresponding function 

f(p,q)=[ fdxdyK(x,y)e- i P^S(q-^-) (7) 

such that (x | / | y) = K(x,y); Weyl ordering has been used[7]. The inverse formula reads 

K(x, y)=J J ^/(p, q y^5(q ^±1) (8) 

and becomes a kind of Fourier transform after integration over q (in order to obtain these 
two formulae one reorders p, q in the exponent in the rhs of (5) using the Baker-Campbell- 
Hausdorff formula and calculates the matrix elements directly). 

The following treatment beyond the WM correspondence also exists: one can avoid ex- 
plicit mentioning the WM correspondence at all and omit the Hilbert space of wavefunctions 
by using the integral transform (7) immediately, then the ^-product of two functions / * g 
corresponds to the convolution of their kernels (i.e. integral transforms(8)): ^ 

K fi , B (x,y) = J d€K f (x,£)K g (€,y). (10) 

In the large noncommutativity limit one can neglect the kinetic term and obtain the 
following equation of motion: ^ = 0, where V ((f)) is considered a polynomial or an analytic 
function (note that the constant term in the lhs is absent). Projection operators are the 
basic thing to solve such an equation. A localized solution will be referred to as a soliton. 

In the next section the general form of any solution radially symmetric w.r.t. spatial 
coordinates is derived. The third section deals with some distribution-valued solutions and 
related topics. The fourth one is reasoning, mainly perturbative and numeric, how does the 
kinetic term affect the physics at finite 6. 



2 General radially symmetric solution 

In[5] radially symmetric projection operators diagonal in the holomorphic representation and 
corresponding functions on the phase space have been worked out. The reason for the present 

Here the same as in (1),(3). In our considerations h — 1. 
2 To verify this one writes the 5-function in (7) in the momentum representation 

S ( q - x -±y)=J^-^ (9) 

or, equivalently, uses (3). 
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examination is that in general adding some non-radially symmetric functions (| m)(n |'s in 
the operator formalism, where m ^ n) or composing a series from them one can in principle 
obtain a radially symmetric result (for example, p 2 + q 2 — r 2 ). In this section we are to 
construct the general ansatz for any radially symmetric solution. In order to do this we 
examine how <j) changes under an infinitesimal coordinate rotation 5q = ep, Sp = —eq: 



(7) can be rewritten as 



5<t>(p,q) = e(p^-q^g) + o(e). (11) 



4>(p,q) = 2 f d£K(q-Z,q + 0e 2lpi . (12) 



After use of (11) one calculates (up to the first order in e) 

+00 

0= f dt{p{d x +d y )K(x,y)-2iqZK(x,y))\ x=q ^ y=q+i e 2lP Z. (13) 

— OO 

Being multiplied by e~ 2ip ^ and integrated w.r.t. p, it becomes 

+00 

0= f dZ(^(d x + d v )K(x,y)d e -2iqZK(x,y))\ x=q - S , v=g+£ 6(Z-e), (14) 



what simplifies to | 



(-\d xx + \x 2 )K{x,y) = (~d yy + l -y 2 )K{x,y). (Hi) 



Separating variables K(x,y) = ^> x (x)'i>y(y) one finds 

{-\d xx + ^x 2 )^ x {x) = EV x (x), 
(-^d yv + ^y 2 )%(y)=E^ y (y). 



(17) 



We are looking for such kernels that convolution of K with K exists (i.e. the integral in 
the rhs of (10) converges). This means that ^ x (and vanishes sufficiently rapidly when 
x — > 00 (correspondingly y — > 00). Similar asymptotic conditions follow from the finiteness 
of energy. ^] There exist some solutions for (17) only when E = (n + -|) and those are the 

wave functions of the Id harmonic oscillator, i.e. the basis of the holomorphic representation 

00 

discussed in [5]. Thus any radially symmetric solution can be represented as C n 4> n , where 

n=0 

C n are some coefficients and <p n are solutions built in[5]. 

3 On some generalized solutions 

Explicit construction. Now we proceed to apply these considerations to the case of 
the SUSY even 4> 4 potential 

V{4>) = -\^ + \<t> i (18) 

The "improper" sign to the quadratic term can describe tachyons as their squared mass is 
negative. On the other hand one can make a constant shift and substitute tf> — ► (j) + const 
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This equation is Fourier self-dual and remains unchanged under the transform 



K(x,y)= j ^K(ky {k * x - k y y \ (15) 



'Thus the WM correspondence is valid. 
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thus expanding the potential near the minimum. In the large noncommutativity limit the 
equation of motion becomes 



-</> + 3 = O. (19) 

It has been pointed out [8] that <f>(x) = ttS 2 (x) satisfies such an equation and[5,8] 

n5 2 (x)*ir5 2 (x) = 1. (20) 

Speaking commonly, how should one understand such a relation? Here we have ir5 2 (x), a 
singular distribution with a compact measure support which can be dealt with as a weak 
limit of a function $(x, e): for instance, 

$(x,e) = ^-e-£,e^+0 (21) 

(momentum representation with a rcgularization factor equal to e~ 2 ). Then carrying 
out the ★-multiplication in the momentum space, one finds that divergences (products of 
singular distributions and their derivatives) cancel each other (what may look like a miracle 
and happens due to the nonlocality of the ^-product). Considering a function like <p(x) — 
wp(x)5(x), where p(x) is some polynomial, one finds that the Fourier transform mapping into 
the momentum space (and, consequently, the ★-multiplication) does not distinguish between 
n p(x) n and n p(0)" . This implies that any function 

0(x) = 7rp(x)S 2 (x) (22) 

such that p(x) is a polynomial and 

p(0) = 1 (23) 

solves (19). At first glance this result seems both evident and trouble-looking, so a little 
discussion on nonlocal product of distributions seems to be necessary. 
First we note that 

d k x (f(x)8{x)) = f(0)d k J(x) (24) 

in a weak sense (as a distribution), where / is a regular function. As an evidence for this one 
can integrate both parts of the above equation with an arbitrary test function and perform 
integration by parts k times. This observation is to prove the dependence only on /(0). 

Another question is to prove the convergence and to calculate the value of the integral 
in the momentum space. No uniform convergence w.r.t. the rcgularization parameter (like e 
in (21)) seems to be valid, but one imposes another regularization in the momentum space 
and successfully calculates the result of the ★-multiplication. This approach is the basic one 
in quantum field theory and the distinction in our case is that no renormalization is needed. 
That is why we say "the divergences cancel each other". So being rigorous, the result in the 
momentum space also converges as a distribution. 

The third thing is the series 

00 

7r5 2 (x) = J2(-l) n <t>n(x) (25) 

ra=0 

converges only as a distribution. For example, replacing (—1)™ with (— e) n , one gets a suitable 
regularization. In particular the vacuum energy becomes finite and equal to 

1 4- 2e 2 

E(e) = - , ' , (26) 
W 4(1 -e 4 )' ' 

but E — ► —00 as £ ^ 1 — 0. Evidence that the configuration C n = (—1)™ is situated at 
the boundary point of the strong convergence range enables one to put forward the equiv- 
alent mathematical approach to the polynomial p(x) as a factor changing the form of the 
regularization. 

The above considerations can also be applied to power series with reasonable asymptotic 
behaviour of the coefficients. 

It is interesting to see what is the corresponding operator in the Hilbert space. From (8) 
one easily calculates K(x,y) — 5(x + y), i.e. this is the spatial reflection P in imaginary 
particle's configuration space. That is the obvious reason why it is no projector. Correspond- 
ing coefficients C n are simply the matrix elements of P in the holomorphic representation: 
C n = (—1)" as the harmonic oscillator wavefunctions have definite parity. For the unity 
operator, of course, C n = 1 and we can construct a projector onto even states 5 (J + P), 
4>=\(l + Tip{x)5 2 (x)). 
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Classical stability. Let us examine whether our solution C„ = (—1)™, = oo is stable. 
For this purpose we expand the energy functional near this configuration up to the second 
order terms and check if the quadratic form to appear is positively defined. No radial 
symmetry of the infinitesimal field change 8<p is implied. The second variation proves to be 
(in operator formalism) |^] 

SE[cf>] = 2 -^-Tr n {-\{Uf + -U^U 2 + 2(^) 2 )), (27) 
5 2 4 

what by use of <j> 2 — I reduces to 

SE[(t>] = ^Tr H (S4> 2 + (<M0) 2 ) > 0. (28) 

9 

Thus the solution is stable. The same considerations apply to <j> = 1, C n = 1. It is clear that 
any solution of (19) to satisfy 

(j)*4>=l (29) 

is classically stable while the solution for which exists n such that C n — is unstable, for 
instance, to the radial fluctuations of the n-th mode because the quartic term gives no rise 
to the quadratic part of expansion now. 

Physical interpretation. First, "localization" of solution is understood comparatively 
as the energy of the <5-like solution is infinite, but it is perfectly localized spatially. 

Coefficients of the polynomial p(x) are unlikely to be important on the classical level and 
seem to be some internal degrees of freedom. The question concerns their meaning. 

In general there exist some different treatments with noncommutative solitons. 

1. Such a localized 2d configuration can be considered to be a D23-brane within D25 
brane[8,9]. The correct brane tension is recovered from the action calculated on the noncom- 
mutative soliton. P{x) can describe excitations living on the world volume of this D23-brane. 

2. In[8] it was proposed that noncommutative solitons are the classical vacua for the 
compact scalar field theory. In our case this statement is supported by the fact that 

E[(j)}\^s 3 {x) = #MU=i = -oo, (30) 
while the energy density remains constant and equal to 

Qe = -^- (31) 

The energy density does not depend on the sign of C n if C n — ±1 at all as <j> n are ★-orthogonal. 
Then such classical vacua can be parametrized by the infinite number of coefficients C n = 
±1. The coefficients of the polynomial p(x) are some excitations. Similar considerations 
concerning the energy density partially apply to the next item. 

Following these ideas, it is interesting to exhibit an example of an instanton interpolating 
between such vacua. To do this we first consider the scalar self-interacting field theory in 

2+1-dimensional Minkowski space-time with the (H ) signature and commutative timclikc 

direction. The action functional is [] 

S[4>] = ^J d 3 x(^/* W ~ V{ff>)), (32) 

from what one reads the energy functional (4) and (now in rescaled noncommuting space 
coordinates) the non-static equation of motion 

W = e A *-^- (33) 

When 8 — oo the first term in the rhs can be neglected. In the euclidean time (or, equiva- 
lently, inverted potential) 

*+ = ™. (34) 

dt 2 86 v ; 



5 We have performed some cyclic permutations under the trace sign. Field variations are considered to be 
small. 

6 No rescaling has been implied yet. 
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Now it is easy to construct an example of an instanton tunneling between two noncommuta- 
tive vacua. Should one restrict the solution to be diagonal at any time moment — oo < t < oo, 
the equation of motion becomes the infinite set of decoupled equations 

^ Cn = C n -Cl (35) 

each of them describes the well known quantum mechanical double well kink if the coefficient 
C n is tunneling from -1 to 1, antikink (1 — > —1) or a constant (1 — > 1 or —1 — > — 1, no 
tunneling) . 

This construction can be applied to the dynamics of radially symmetric solitons and the 
equations of motion are written in the Minkowski space-time. Recently some new results 
concerning noncommutative soliton scattering[10] were obtained using the adiabatic approx- 
imation. [] The typical feature is that the equations (35) are identical for all modes. If all 
the initial conditions are the same (up to a sign), the solution is (f>(r,t) — C(t)<f>(r), i.e. the 
time dynamics reduces to the scaling factor in front of the soliton. As V(4>) — > oo when 
4> — > oo, all the coefficients C n remain finite. 

If small excitations mentioned above are radially symmetric, they can be described as 
small oscillations of C n near the minimum of potential. 

3. Some vortex- like solutions in noncommutative abelian Higgs model (which has similar 
potential for the scalar field) have recently been discussed[ll]. In this context (and also 
for the brane interpretation) the subtlety to arise is the uniform energy density while it is 
significantly singular for a commutative brane or a vortex. It is not completely clear for what 
vortex's internal degrees of freedom p(x) is responsible. 



4 Behaviour of solutions with finite 9 

Recurrent relation. The question still unanswered completely is how does the kinetic 
term affect the physics at finite 6. To examine some aspects one should write the equations 
of motion at finite 6. It is known[5,I2] that in the operator formalism 

i^d^M]^], (36) 

where a\ a are creation and annihilation operators defined in Appendix. Performing varia- 
tion in (4) w.r.t. </>, one obtains 

3V - 

[a\[a,<j>]] + U,a%a] = -e— {<!>). (37) 

By use of Jacobi identity and commutation relations for creation and annihilation operators 
lhs yields 

[a\[a^]] = -\e^). (38) 
With the radially symmetric ansatz the equation of motion simplifies to the recurrent relation 

a i , 

(n + l)C„+i - (2n + 1)C„ + nC n -i = -0V (C n ), n > 0, 

(39) 

c 1 -Co = -ev\c ),n = o. 

Should one try to find an appropriate potential V(4>) for a solution C n — (—l) n A to exist, 
it will end in no positive result as in this case (39) reduces to 

2(-I) n+1 (2n + 1)A = -V'((-l) n A). (40) 

As it has been emphasized, an equation like (39) is hard to solve, but some effects could be 
studied qualitatively or numerically. 



7 The calculations below are closed for 8 — oo. 
8 For derivation see Appendix. 
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Qualitative effects. l.The differential equation for this scheme is 



nC" +C' = ^9C(C 2 - 1). (41) 

A similar one occurs in the construction of the usual commutative Abrikosov-Nielsen-Olesen 
vortex and the only acceptable asymptotic for us is that the coefficients C (n) = C n vanish at 
the infinity. The mode expansion (j> = ^2 C n cf) n is somewhat like a discrete Fourier transform. 
Then the singular (or at least localized) vortex-like solution in the discrete momentum space 
should lead to derealization in the field 0-space. Thus, the (5-like solution becomes finite-size. 
The effect of finite 9 is imposing a rcgularization. 
2. The energy functional is[12] 

E W = \ £ « 2n + l ) C l - 2(n + l)C n+1 C n + 6V(C n )). (42) 

^ n=0 

Analyzing the first two terms as a perturbation, i.e. simply calculating the energy functional 
for the solution at 9 = oo one observes that the second addend in the kinetic term causes 
the level splitting and breaks the 2°°-times degeneration. Perfectly localized <5-like solution 
becomes the highest energy level because all the perturbative terms in (41) appear to be 
positive for the coefficients with alternating sign. This is the reason to understand of what 
importance is the (5-like radially symmetric solution. To do that one observes that the results 
for 9 = oo are true for any set of mutually orthogonal projectors. It is the kinetic term that 
causes the [/(oo) symmetry breakdown. For an arbitrary operator like 

= UUiag({X n })U, (43) 
with U being unitary the energy functional is equal[12] up to the coupling constant 

oo oo oo oo 

£[{A„},{lW}] = ]TA2(l + 2^m|[/ m „| 2 )-2 ^mK\A mn \ 2 + 9Y J V(\n), 

n— m— m.n—0 n—0 

(44) 

where 

oo 

A mn = ]Tv / k/ fc „t4*_ ljm , (45) 
fc=i 

If we hold {A„} fixed, then the first order variation of E w.r.t. U is 

„ dE dE - 

oE = — — \u=iOU mn + —= — \u=ioU mn = 

(46) 

2^(3n+l)(<JE/ nn + Ji/ nn ) = 

n=0 



as in the first order 



SU nn + 6U nn = (47) 



for U is unitary. 

Thus for 9 = oo the matrix U serves as coordinates along the vacuum valley, but at finite 
9 the special significance of radially symmetric solutions becomes clear. 



Numeric study. Our method is different from that in[12]. It is clear that specifying Co 
defines all the other coefficients and provides an opportunity for the numeric study. In our 
study 9 = 100 and remains constant. The main goal is to see how does a finite-^ solution 
look like. Some graphs are presented below for Co = 10~ 10 . 

The equations are strongly nonlinear, so as we start from Co = 10~ 10 it first grows in ten 
orders up to C <~ 1 and then slowly converges to (see the first graph). In the second graph 
we show the sum of the 100 first terms in the series for 4>(r). The third one illustrates that 
4>(0) remains finite and thus the solution is finite-size. 

Another interesting phenomenon is that a small change in Co can make the solution 
divergent. For example, there exist at least two disconnected ranges of convergence Co < 
3.82 • 10~ 10 and 9.1 • 10~ 10 < C < 11.6 • lO" 10 , 9 = 100. This point agrees well with the 
soliton hierarchy first observed in [12]. 
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Dependence C(n) 




±0.4 



Solution in configuration space: sum of the first 100 terms 




5 Concluding remarks 



In the present paper we consider 2+1-dimensional noncommutative classical scalar field the- 
ory. The standard and basic tool for solving equations of motion at 9 — oo is the Weyl-Moyal 
correspondence between the product of operators acting on the Hilbert space and the nonlo- 
cal product of functions. We mention how one can solve the equations of motion in principle 
without referring to the Weyl-Moyal correspondence. Then the general ansatz for a radially 
symmetric solution is obtained confirming the celebrated results of [5] . Next these results are 
applied to studying generalized (distribution- valued) solutions and their physical assignment. 
The ^-product of distributions is defined as a distribution in the momentum space. Small 
excitations are shown to change the regularization. An approach to the accurate study of the 
radially symmetric solution dynamics is given and the effect of the time dynamics coming 
to the field rescaling is described. Then we proceed to examine the physics at finite 9. The 
effects of delocalization and taking down the degeneration are analyzed qualitatively and 
numerically. The special role of radially symmetric solutions is clarified. 

There arise different questions of interest such as dynamics of fields without radial sym- 
metry and general scattering problem. It is interesting to take quantum corrections into 
account. Among other problems there are coupling with the gauge field and searching out 
some closed form solution and analysis for finite 9. 
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7 Appendix 

Derivation of the recurrent relations. In the holomorphic representation the key 
role is played by the two operators 

-t q + w 
a = — 

Commutation relations for them [a, = I do have a representation in the space of analytic 
functions f(z) with the inner product[13] 

(A,/2>= / ^e(-*>A(i)/ a (*). (49) 



2tti 

The orthonormal state vectors are represented as 



n) i> n (z) = ^4, (50) 



creation and annihilation operators as 



<V <-> z. 



(51) 



a «-> d. 

An arbitrary operator U with matrix elements equal to 

U mn = {tpm\U\ $n) (52) 
can be represented as the integral operator with the kernel f\ 

A{z,z) =y2u mn -=-=, (53) 



3 A(z,w) is an analytic function (some series) of two complex arguments, and in general z ^ w. 
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(Uf)(z) = J ^-e(-^A(z,z)f(z). (54) 

Now all the needed commutators in the equation of motion can be simply computed using 
the formulae 

a^U «-> zA(z,z), 
aU <-> dA(z, z), 

. , (55) 
C/a T «-> dA(z, z), 

Ua «-> zA(z, z). 

Under the assumption of the radial symmetry C/ mn = C m 5 mn the result is 

- zd - zB + zz - l)A(z, z) = -^ -gr(C n )i) n {z)ip n (z), (56) 

n— 

and after the simple direct calculation of matrix elements one arrives to (39). The energy 
functional is obtained in a similar way. 
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